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Abstract 

This review describes the main experimental facts and a number of theoretical models concerning 
the pseudogap state in high - temperature superconductors. On the phase diagram of HTSC - 
cuprates the pseudogap state is observed in the region of current carrier concentrations less than 
optimal and leads to a number of anomalies of electronic properties. These pseudogap anomalies 
apparently are due to fluctuations of antiferromagnetic short - range order developing as system 
moves towards antiferromagnetic region on the phase diagram. Electron interaction with these 
fluctuations leads to anisotropic renormalization of electronic spectrum and non - Fermi liquid 
behavior on certain parts of the Fermi surface. We also discuss some simple theoretical models 
describing the basic properties of the pseudogap state, including the anomalies of superconducting 
state due to pseudogap renormalization of electronic spectrum. 
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INTRODUCTION. 



The studies of high - temperature superconductivity (HTSC) in copper oxides remains 
one of the central topics of the condensed - matter physics. Despite of the great effort of both 
experimentalist and theorists the nature of this phenomenon is still not fully understood. 
It is common belief that the basic difficulties here are due to quite unusual properties of 
these systems in the normal (non - superconducting) state, while without the explanation of 
these anomalies it is hard to expect the complete understanding of microscopic mechanism 
of HTSC. In recent years one of the main directions in physics of HTSC were the studies 
of the so called pseudogap state M, which is observed in the region of the phase diagram, 
corresponding to current carrier concentrations less than the optimal concentration (i.e. that 
corresponding to the maximal critical temperature of superconducting transition T c ). This is 
usually called the "underdoped" region. In this region a number of anomalies of electronic 
properties is observed both in normal and superconducting states, which are interpreted 
as due to a drop of the density of single particle states and anisotropic renormalization 
of the spectral density of current carriers. Understanding of the nature and properties of 
the pseudogap state remains the central problem in any approach towards explaining the 
complicated phase diagram of HTSC - systems. This problem was already discussed in some 
hundreds of experimental and theoretical papers 1 . 

The aim of this review is to report basic experimental facts on the observation of the 
pseudogap state in underdoped HTSC - cuprates, as well as the discussion of a number of 
simple theoretical models of this state. The review is not thorough, both in describing the 
experimental data and also in theoretical part. Experiments are discussed rather briefly, 
taking into account the existence several good reviews |fj, ||, [|, || |[ . Theoretical discussion 
is also rather incomplete and reflects mainly the personal position of the author. There exist 
two main theoretical scenarios for the explanation of pseudogap anomalies in HTSC - sys- 
tems. The first is based upon the model of Cooper pairs formation already above the critical 
temperature of superconducting transition |7|, ||, |], 03 , with phase coherence appearing for 

1 Note that only the cond-mat electronic preprint archive contains more than 600 papers somehow related 
to the physics of the pseudogap state, while at the last major M 2 S — HTSC — VI Superconductivity 
Conference in Houston (February 2000) four oral sessions were devoted to this problem - more than to 
any other problem in the field of HTSC. 
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T < T c . The second one assumes that the appearance of the pseudogap state is due to 
fluctuations of some short - range order of "dielectric" type which exist in the underdoped 
region of the phase diagram. Most popular here is the picture of antiferromagnetic (AFM) 
fluctuations JTTJ, |T2|, [TB|, [T4], TE\ , while the similar role of fluctuations of charge density waves 
(CDW), structure or phase separation on microscopic scale can not be excluded. The author 
believes that some recent experiments strongly suggest that this second scenario is correct. 
Accordingly, in theoretical part of this review we limit ourselves only to discussion of this 
kind of models with detailed enough explanation of the results obtained by the author and 
his collaborators, as well as some other authors whose works are somehow close to ours in 
main ideas and in approaches used. Thus, we mainly assume the relevance of the model of 
antiferromagnetic fluctuations, though it is still can not be considered as fully confirmed. 

The references are also apparently incomplete, moreover these do not reflect any con- 
siderations of priority. It is assumed that the reader will find the additional references in 
papers cite here. The author expresses his excuses in advance to those numerous authors, 
whose papers are not cited here, mainly not to make the reference list too long. 



BASIC EXPERIMENTAL FACTS. 

Typical variants of the phase diagram of HTSC - cuprates are shown in Fig.[l|. Depending 
on concentration of current carriers (in most cases holes) in highly conducting C11O2 plane 
a number of phases is observed as well as regions with anomalous physical properties. For 
small concentration of holes all known HTSC - cuprates are antiferromagnetic insulators. As 
hole concentration grows Neel temperature Tjy drops rapidly from the values of the order of 
hundreds of K, becoming zero at hole concentration p of the order of 0.05 or even less, while 
the system becomes metallic (though this metal is rather bad). With the further growth 
of hole concentration the system becomes superconducting and superconducting critical 
temperature at first grows with hole concentration, becoming maximal at po ~ 0.15 — 0.17 
(optimal doping), then drops to zero at p ~ 0.25 — 0.30, while in this (overdoped) region 
the metallic behavior persists. Metallic properties in the region of p > p^ are more or less 
traditional (Fermi - liquid), while for p < po the system is a kind of anomalous metal, not 
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FIG. 1: Two possible types of the phase diagram of HTSC - cuprates. 

described (in the opinion of majority researchers) by Fermi - liquid theory 2 . 

Anomalies of physical properties which are attributed to the appearance of the pseudo- 
gap state are observed in metallic phase for p < po and temperatures T < T*, where T* 
falls from temperatures of the order of Tjy for p ~ 0.05, becoming zero at some "critical" 
concentration p c , which is may be slightly larger than p (Fig.|](a)), e.g. according to Ref. 
this happens at p = p c « 0.19. In the opinion of a number of authors (mainly supporters 
of superconducting scenario of pseudogap formation) T* gradually approaches the line of 
superconducting T c close to the optimal concentration p (Fig.[L](b)). Below we shall see that 
the majority of new experimental data apparently confirm the type of the phase diagram 
shown in Fig.|l](a) (see details in Ref. 0). It should be stressed that temperature T*, in 
the opinion of majority of researchers, does not have the meaning of a temperature of some 
kind of phase transition. It only determines some characteristic temperature scale below 
which the pseudogap anomalies appear in different measurements. Any kind of singularities 
of thermodynamic characteristics in this region of phase diagram are just absent 3 . The 

2 Note that the question of presence or absence of Fermi - liquid behavior in HTSC - systems is rather 
obscured in discussions of numerous authors standing on opposite positions and even using different 
definitions of Fermi - liquid. In the following we shall use this terminology according to the opinion of 
rather badly defined "majority" . 

3 Note, however, an opposite opinion expressed in a recent paper where the line of T* on the phase 




hole concenlrelion, p 
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general statement is that all these anomalies can be related to some drop (in this region) 
in the density of single - particle excitations close to the Fermi level, which corresponds to 
a general concept of the pseudogap 4 . In this case the value of T* is just proportional to 
energy width of the pseudogap. Sometimes an additional temperature scale T 2 * is also intro- 
duced, as shown in Fig. [TJ(b), attributed to a crossover from "weak" to "strong" pseudogap 



regime [H|, usually referring to some change in spin system response in the vicinity of this 
temperature. In the following we shall not be dealing with these details. 

Now let us describe the most typical experimental signatures of pseudogap anomalies in 
HTSC - cuprates. 

Specific Heat and Tunneling. 

Consider first the data on electronic specific heat of HTSC - cuprates. In metals this 
contribution is usually represented as C = j(T)T, so that in the normal state (T > T c ) 
7 ~ N(0), where N(0) - is the density of states at the Fermi level. For T = T c we have 
the well known anomaly due the second order transition and the value of j(T) acquires 
the characteristic peak ("jump" or discontinuity). As an example in Fig.|2| we show some 
typical data obtained for YQ & CaQ,2Ba2Cu^O-js with different values of 5 [18]. In optimally 
doped and overdoped samples j(T) remains practically constant for all T > T c , while for 
underdoped samples a significant drop of 7(T) is observed for temperatures T < 150 — 200K. 
This observation directly signifies to the drop in the electronic density of states at the Fermi 
level and formation of the pseudogap for T < T*. 

Note also that the value of the "jump" of specific heat at T c is much reduced as we go 
into the underdoped region. More detailed analysis |J shows that the appropriate value of 
A'jc drops rather sharply, starting from the "critical" concentration of carriers p c ~ 0.19, 
which is considered to be the point of appearance of the pseudogap in electronic spectrum. 



diagram was directly connected with some "hidden" symmetry breaking, strongly smeared by internal 
disorder. 

4 The concept of pseudogap was first formulated by Mott in his qualitative theory of disordered (non - 



crystalline) semiconductors 17 1 . According to Mott the term "pseudogap" describes a region of diminished 
density of states within the energy interval corresponding to a band gap in an ideal crystal, representing 
the "reminiscence" of this band gap which is conserved after strong disordering (amorphization, melting 
etc.) 
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FIG. 2: Electronic specific heat coefficient 7 in overdoped (a) and underdoped (b) 
Y . 8 Ca . 2 Ba 2 Cu 3 O 7 _s @- 

Using more or less traditional expressions of BCS theory, specific heat data can be used to 
estimate the value of superconducting energy gap in electronic spectrum and its temperature 
dependence. In Fig.|3] we show the data of Ref . |1| , obtained for Y Ba 2 Cu^Ojs (YBCO). 



It is seen that thus defined energy gap does not go to zero for T = T c (as it should be 
in a traditional superconductor) and there is a significant "tail" in the region of higher 
temperatures. This effect is stronger for more underdoped samples. Often these data are 
naively interpreted as signifying Cooper pair formation at temperatures T > T c . 

Pseudogap formation is clearly seen in the experiments on single particle tunneling. 
In a widely cited Ref.[|20| tunneling measurements were performed on single crystals of 
Bi 2 Sr 2 C 'aCu 2 Os+s (Bi — 2212) with different oxygen content. For underdoped samples the 
pseudogap formation in the density of states was clearly observed for temperatures higher 
than T c . This pseudogap smoothly evolved into superconducting gap for T < T c , which 
again was seen by many as a direct confirmation of superconducting nature of the pseudo- 
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FIG. 3: Temperature dependence of the energy gap for samples of Y Ba<iCuj,0§+ x with different 
compositions, determined from specific heat measurements |l9[| . The number at different curves 
denote the appropriate values of oxygen content x. 

gap. Some signs of pseudogap existence was observed in this work also for slightly overdoped 
samples. 

While discussing tunneling experiments in HTSC - cuprates, there is always a question 
of the surface quality of the samples. Thus, especially interesting the recent studies of Refs. 

|22| , where an intrinsic tunneling was measured on the so called mesa - structures 5 , 
specially formed (by microlitography) on the surface of the same system Bi — 2212. These 
experiments clearly demonstrated the existence of superconducting gap approaching zero at 
T = T c on the background of smooth pseudogap, which exists at higher temperatures. The 
appropriate tunneling date are shown in Fig.f|. Typical peaks due to superconducting gap 
opening are seen on the background of smooth minimum in the density of states due to pseu- 
dogap. It is especially important that in Ref.|22[ it was demonstrated that superconducting 
peculiarities of tunneling characteristics were destroyed by the external magnetic field, while 
the pseudogap remained practically independent of the field strength, which clearly shows 



5 "mesa (Sp.), a land formation having a flat top and steep rock walls: common in arid and semiarid parts 
of the U.S. and Mexico" (The Random House Dictionary of the English Language. Random House, Inc., 
1969.). 
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FIG. 4: Differential conductance of nearly optimally doped Bi — 2212 for different temperatures. 



At the insert - results for higher T [21]. 



its non - superconducting nature. Appropriate data are shown in Fig| 



NMR and Transport Properties. 

Formation of the pseudogap signals itself also in transport properties of HTSC - systems 
in normal state, Knight shift and NMR - relaxation. In particular it is thought to be a 
reason for the change of the standard (for optimally doped samples) linear temperature 
dependence of resistivity in the region of T < T* for underdoped samples. The value of 
the Knight shift in underdoped samples becomes temperature dependent - for T < T* it 
drops rather rapidly. Analogous behavior is observed in underdoped region for the value of 
(TTx) -1 , where 7\ - is NMR relaxation time. Let us remind that in usual metals the Knight 
shift is proportional to N(0), while (TTi) -1 ~ A^ 2 (0) (Korringa relation). Resistivity p is 
usually proportional to the scattering rate (inverse mean - free time) 7 ~ N(0). Thus, the 
notable suppression of these characteristics is naturally explained by the drop in the density 
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FIG. 5: Differential conductance of nearly optimally doped Bi — 2212 for different temperatures 
in magnetic field H = (thin lines) and for H = 14T (thick lines) [p2j]. 

of states N(0) at the Fermi level 6 . Note that this interpretation is obviously oversimplified, 
moreover when we discuss the temperature dependence of these physical characteristics. For 
example, in the case of resistivity this dependence is usually mainly determined by inelastic 
scattering processes and the physical mechanism of these in HTSC - cuprates is still far 
from being clear. The drop in the density of states, interpreted as partial dielectrization of 
electronic spectrum, may also lead to the appropriate growth of resistivity. 

In Fig.|6] taken from Ref. M we show the summary of experimental data for under- 
doped samples of Y Ba^Cu^O^ using the results of Refs. |23|, p4] , p5fl . In Fig.f7| we show 
the data of Ref. [26] on the measurements of (TTi) -1 on 63 Cii-nuclei in slightly underdoped 



YBa2Cu 3 7 _s m strong enough magnetic field. It is seen that the magnetic field dependence 
of the data is practically absent, giving the strong support for non - superconducting nature 



6 Historically the first data on the significant drop of the density of states in underdoped cuprates were 
obtained in both NMR and magnetic neutron scattering. Accordingly at this initial stage the term 
"spin gap" was introduced and actively used. However, from further studies it became clear that the 
analogous effects are observed not only for spin degrees of freedom and the term "pseudogap" became 
widely accepted. 
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FIG. 6: Influence of the pseudogap formation on resistivity, Knight shift and NMR relaxation in 
underdoped Y Ba2Cu/\Q%. The changes in temperature dependences in the region of T < 200 — 
300i^ is attributed to pseudogap formation in the density of states M. 
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FIG. 7: NMR relaxation rate (TTi) -1 for slightly underdoped Y Ba2Cu^Ojs i n strong magnetic 
field, directed along the c-axis of the sample [^]. Vertical line shows T c in zero field . The straight 
line parallel to the data estimates the expected shift of data in the case of pseudogap suppression 
by magnetic field, its position is determined by the observed shift of T c . 

of the pseudogap, as in the opposite case magnetic field should significantly change the size 
of the effect. Note that the expected influence of magnetic field upon the Knight shift and 
NMR relaxation of m Cu was observed in slightly overdoped TlSr 2 CaCu20e^ where it was 
accurately described by the effects of suppression of superconducting fluctuations |27j . This 
fact can be seen as an evidence of disappearance of the pseudogap of non - superconducting 
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FIG. 8: Dependence of energy scale characterizing the width of the pseudogap E g in YBCO on 
concentration of holes following from different experiments ||. 

nature (independent of magnetic field) at some lower concentration of current carriers. 

Let us stress that in different experiments the characteristic temperature T* below which 
the anomalies attributed to pseudogap appear can somehow change depending on the quan- 
tity being measured. However, in most cases there is a systematic dependence of T* on 
doping and this temperature tends to zero close to optimal concentration of current carriers 
or at slightly higher concentrations. It is clear that T* does not represent any strictly de- 
fined temperature (e.g. of some phase transition) but only determines some characteristic 
crossover temperature (energy) scale below which pseudogap anomalies appear . This en- 
ergy scale may be determined from experimental data using some rough model of the energy 
dependence of the density of states close to the Fermi level within the pseudogap region (e.g. 
triangular with the width of E g , giving T* = 0AE g ) ||. In Fig.|8| we present the summary of 
the values of pseudogap width determined in this way from different experiments on YBCO 
for different hole concentrations 0]. It is seen that the pseudogap "closes" at some critical 
value of p c ~ 0.19. Once again let us stress that both the line of T* on the phase diagram 
and "critical" concentration p c apparently possess only qualitative meaning. However, there 
are also attempts to interpret p c as some kind of a "quantum" critical point |2S[. 
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FIG. 9: Real part of optical conductivity in ab - plane of Lai^SrQ^CuO^ for different tempera- 
tures [3C]. 

Optical Conductivity. 

Pseudogap formation in underdoped HTSC - cuprates is also clearly seen in numerous 
measurements of optical conductivity, both for electric field vector polarized along the highly 
conducting C11O2 plane and along orthogonal direction of tetragonal axis c. These experi- 
ments are reviewed in detail in Ref. pi f2^]. As a typical example in Fig.|9| we show the data 
of Ref. |30t1 on optical conductivity in the CuOi - plane of underdoped Lai^Sr^^CuO^. 



Characteristic feature of these data is the presence of rather narrow "Drude - like" peak 
in the frequency region u < 250cm" 1 and the appearance of "pseudogap suppression" of 
conductivity in the interval of 250 — 700cm" 1 accompanied by smooth maximum around 
uj ~ 800cm" 1 . Especially sharply pseudogap anomalies are seen after the experimental data 
are expressed in terms of the so called generalized Drude formula || 29fl , which introduces 



the effective (frequency dependent) scattering rate and effective mass of current carriers: 
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FIG. 10: Generalized scattering rate and effective mass as functions of the frequency of the external 



field obtained from the fit of optical data for Lai^Sro^CuO^ to generalized Drude formula [30 1: 
(a), (c) - for temperatures T < 250K, when the linear growth of 1/r with frequency is observed 
for uj > 700cm _1 , (b), (d) - higher temperatures, when the linear growth of 1/r is practically 
invisible. 

where a{uS) - the observed complex conductivity, u p - plasma frequency, m - mass of the 
free electron. Note that the use of expressions fll]), (0) is just another form of expression 
of the same experimental data in the form of 1/r (a;) and m*{uo) instead of two more usual 
characteristics — the real and imaginary parts of conductivity. Accordingly there is no 
especially deep physical meaning in both t(oj) and m*(uj). However, this form of presentation 
of experimental data in optics is widely used in the literature. 

In Fig.[Tl] we present the results of such a representation of experimental data of the sort 
shown in Fig.| @. It is seen that the effective scattering rate for temperatures less than 
T* ~ A50K is sharply suppressed for frequencies of an external field less than 700cm _1 , 
while for higher frequencies grows linearly with u, demonstrating the anomalous non 
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- Fermi liquid behavior. The presence of this dip in 1/t(u) is most vivid reflection of the 
pseudogap in optical data. At the same time we have to note that the frequency dependence 
of m*(u) shown in Fig.[10] apparently has no clear interpretation. 

Analogous behavior of 1/t(u) is also observed in data on transverse optical conductivity 
a c (uj), practically for all HTSC - cuprates studied in underdoped region |2|, |2Pj . 



Fermi Surface and ARPES. 



Especially spectacular effects due to formation of the pseudogap are observed in exper- 
iments on angle - resolved photoemission (ARPES) 0, j5|. This method is, first of all, 



decisive in studies of the topology of the Fermi surface of HTSC - cuprates j|, [31], [32]], being 
practically the sole source of this information. 

In Fig.[ll| we show the data of Ref.Q] on the Fermi surface of La2- x Sr x CuO± for two 
different compositions. In overdoped state (x = 0.3) Fermi surface is electronic and centered 
around the point T(0,0) in Brillouin zone. As concentration x diminishes the topology 
of the Fermi surface changes, so that for optimal concentration of Sr and in underdoped 
region (x = 0.1) it becomes hole - like and centered around the point Y(tt, n). Precisely this 
topology is observed in majority of ARPES - experiments for almost all HTSC - cuprates. 



As an interesting example in Fig.[12] from Ref. |34| we show the hole - like Fermi surface of 
most extensively studied (by ARPES) system Bi — 2212 in overdoped state. The spectacular 
feature here is the observation of large flat parts of the Fermi surface, orthogonal to direction 



YM and symmetric directions. Independently this result was confirmed in Ref.|55| for 
optimally doped system. The presence of the flat parts (patches) on the Fermi surface may 
be very important for the construction of microscopic theories of electronic properties of 
HTSC - systems. Dropping these fine details in the first rough enough approximation the 
observed topology of the Fermi surface and the spectrum of elementary excitations in the 
CuOi plane can be rather accurately described by the tight - binding model: 

£k = — 2t(cos k x a + cos k y a) — At' cos k x a cos k y a (3) 

where t ~ 0.25eV - transfer integral between nearest neighbors, and i! - transfer integral 
between second nearest neighbors, which may change between t' ~ — 0A5t for Y 'BCL2CU3O7-S 
and t' — 0.25i for La,2- x Sr x CuO '4, a - lattice constant. 
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FIG. 11: Fermi surface of La,2- x Sr x CuO\ for overdoped (a) and underdoped (b) cases 33] 



Gatt et al., Fig. 3 




FIG. 12: Fermi surface of overdoped Bi — 2212 jj^j]. Circles represent points where measurements 
were performed. At the inserts - energy distribution of photoelectrons along the directions a and 
b, shown by arrows. 
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Quite recently the simple picture described above was seriously revised in Refs. |35], |36], |37|, 
|38| using the new ARPES data obtained with synchrotron radiation of higher energy, than 
in previous experiments. It was claimed that the Fermi surface of Bi — 2212 is electronic 
and centered around point V. The main discrepancies with previous data are in the most 
interesting (for us) vicinity of the point (0, it) in the inverse space. These claims, however, 



were met with sharp objections from other groups of ARPES experimentalists [EM EDI [IT 



42| , [43|]. The question remains open, but in the following we shall adhere to traditional 
interpretation. 

There are several good reviews on the observation of pseudogap anomalies in ARPES 
experiments || [|, ||, so below we shall limit ourselves only to basic qualitative statements. 
The intensity of ARPES (energy distribution of photoelectrons) in fact is determined as [|J : 

/(kw) = J (k)/HA(ku;) (4) 

where k - is the electron momentum in the Brillouin zone, u - the energy of initial state 
with respect to the Fermi level (chemical potential) 7 , Io(k) includes kinematic factors and 
the square of the matrix element of electron - photon interaction (and in rough enough 
approximation can be considered a constant). The value 

A(k, u) = --ImG(k, u + i6) (5) 

7T 

where G(k, u) - is the Green's function, represents the spectral density of electrons. Fermi 
distribution f(ui) = [exp(u/T) + 1] _1 reflects the fact that only electrons from the occu- 
pied states participate in photoemission. Thus in a rough approximation we may say that 
ARPES experiments directly measure the product of f(u)A(kuj), and we can get the direct 
information on the spectral properties of single - particle excitations in our system. 

In FigfTBI we show typical ARPES data for Bi 2 Sr 2 CaCu20s H, measured in three 



different points on the Fermi surface for different temperatures. The appearance of the 
gap (pseudogap) is reflected in the shift (to the left) of the leading edge of photoelectron 
energy distribution curve, compared to the reference spectrum of a good metal {Pt). It is 
seen that the gap is closed at different temperatures for for different values of k, while its 
width diminishes as we move away from (0, 0) — (0, 7r) direction. Pseudogap is completely 



7 In real experiments ui is measured with respect to the Fermi level of some good metal like Pt or Ag which 
is put into electric contact with the sample. 
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FIG. 13: ARPES spectra in three different points a, b, c on the Fermi surface of Bi — 2212 for 
underdoped sample with T c = 85K. Thin lines - ARPES spectrum of a reference sample of Pt 



absent in the direction of the zone diagonal (0, 0) — (tt, ii). At low temperatures these data 
are consistent with the picture of (i-wave pairing, which is well established from numerous 
experiments on HTSC - cuprates ffEfl . It is important, however, that this "gap" in 
ARPES data is observed also for temperatures significantly higher than the temperature of 
superconducting transition T c . 

In Fig.pJ] we show the angular dependence of the gap width in the Brillouin zone and 
the temperature dependence of its maximum value for the number of samples of Bi — 2212 
with different oxygen content, obtained via ARPES in Ref . ||48|| . It is seen that with the 
general <i-wave symmetry, the gap in the spectrum of optimally doped system becomes zero 
practically at T = T c , while for underdoped samples there appear characteristic "tails" 
in the temperature dependence of the gap for T > T c , quite similar to those shown in 
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FIG. 14: Angular (a) and temperature (b) dependence of energy gap in Bi — 2212, obtained from 
ARPES data || for samples with T c = 87 K (nearly optimally doped), T c = 83K and T c = WK 
(underdoped). 

Fig|| Qualitatively we may say that the formation of anisotropic (in reciprocal space) 
pseudogap for T > T c , continuously evolving into superconducting gap for T < T c , leads to 
a "destruction" of the parts of the Fermi surface of underdoped samples already for T < T* 
and around the point (0, 7r) (and symmetric to it), and that the effective width of these parts 
grows as the temperature is lowered |49| . 

Central point is, of course, the evolution of spectral density A(k F , uj) at the Fermi surface. 



Under rather plausible assumptions this may be directly obtained from ARPES data |4J| . In 
case of electron - hole symmetry we have A(kp, uj) = A(kp, —uo) (which is always valid close 
enough to the Fermi level, in real systems for \uj\ less than few tenths of meV), so that, with 
the account of /(— to) = l — f(—uj), from (|) for k = kp we immediately get + I(—u) = 
A(k.F,co). Thus, the spectral density at the Fermi surface may be obtained directly from 
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FIG. 15: Symmetrized ARPES spectra for overdoped sample of Bi — 2212 with T c = 82K (a) 
and underdoped sample with T c = 83K (b) at the point of intersection of the Fermi surface with 
Brillouin zone border (0, 7r) — (vr, tt) 



experiments, more precisely from the symmetrized ARPES - spectrum I(uj) As an 

example of such approach, in Fig.[T5| we present the data of Ref. [TOj for underdoped sample of 
Bi — 2212 with T c = 83K and overdoped sample with T c = 82K for different temperatures. 
It is seen that the existence of the pseudogap is reflected in typical two - peak structure 
of the spectral density, appearing for the underdoped system at temperatures significantly 
higher than T c . 

Let us stress that well defined quasiparticles correspond to sharp enough peak of the 
spectral density A(kp, u) at uj = 0. This type of behavior is practically unobservable in 
HTSC - cuprates, at least for T > T c . In principle, this does not seem surprising - it is 
rather difficult to imagine well defined quasiparticles almost in any system for T > lOO/H 
However, the much improved resolution of modern ARPES installations apparently allows to 
conclude, that the actual width of the peak observed is larger than experimental resolution, 
so that the problem can be studied experimentally plf . It appears that in superconducting 
state, for T T c , in the vicinity of the crossing of the Fermi surface with the diagonal 
of the Brillouin zone ((0,0) — (it, it) direction), there appears sharp enough peak of the 
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spectral density, corresponding to well denned quasiparticles [^1|]. However in the vicinity 
of the point (0, it) the Fermi surface is "destroyed" by superconducting gap, corresponding 
to rf-wave pairing, which leads to the two - peak structure of the spectral density. 

The smooth evolution of ARPES pseudogap observed for T > T c into superconducting 
gap corresponding to d-wave pairing for T < T c is often considered as an evidence of super- 
conducting nature of the pseudogap. However, most probably this is wrong. In this respect 
we should like to note the very interesting paper ||52| , where the ARPES measurements were 



performed on dielectric oxide Ca^CuOiCli-, which is structurally analogous to La 2 CuO^ and 
becomes a high - temperature superconductor after doping by Na or K. ARPES measure- 
ments on this system were possible in insulating phase due to especially good quality of the 
sample surface. The exciting discovery made in Ref. |5l| was apparently first observation of 



the "remnant" Fermi surface in this Mott dielectric, which is everywhere "closed" by the 
gap (apparently of the Mott - Hubbard nature). At the same time a strong anisotropy of 



this gap in inverse space with the symmetry of <i-wave type was also observed [52], very 
similar to analogous data for HTSC - systems in metallic phase. It seems plausible that this 
anisotropic gap has the same nature as high - energy pseudogap say in Bi — 2212 |52|]. Of 
course, it is quite clear that there is no any kind of Cooper pairing in this insulator. 

Other Experiments. 

Pseudogap behavior is also observed in a number of other experiments, such as electronic 



Raman scattering [53] and magnetic scattering of neutrons |54|. All these experiments are 
interpreted as an evidence for significant suppression of the single - particle density of states 
close to the Fermi level of underdoped HTSC - cuprates for temperatures T < T*, that is 
already in normal phase. Due to a space shortage we shall not discuss these data in any 
detail, moreover that a good review of these can be found in papers cited above, as well as 
in Ref. [01. 
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THEORETICAL MODELS OF PSEUDOGAP STATE. 



Scattering by Fluctuations of Short - Range Order: Qualitative Picture. 

As was already noted above there exist two main theoretical scenario to explain pseudogap 
anomalies of HTSC - systems. The first is based on the model of Cooper pair formation 
above the temperature of superconducting transition (precursor pairing), while the second 
assumes that formation of the pseudogap state is due to fluctuations of short - range order of 
"dielectric" type (e.g. antiferromagnetic or charge density wave) existing in the underdoped 
region of the phase diagram and leading to strong scattering of electrons and pseudogap 
renormalization of the spectrum. This second scenario seems more preferable both because 
of experimental evidence 8 , described above, and due to the simple fact that all pseudogap 
anomalies become stronger for stronger underdoping, when system moves from optimal (for 
superconductivity) composition towards insulating (antiferromagnetic) region on the phase 
diagram. 



Consider typical Fermi surface of current carriers in C1LO2 plane, shown in Fig. 16. Phase 



transition to antiferromagnetic state leads to doubling of lattice period and appearance of 



"magnetic" Brillouin zone in reciprocal space, shown by dashed lines in Fig.|16] . If the 
energy spectrum of current carriers is given by (|3]) with t' = 0, then for the case of half - 
filled band the Fermi surface is just a square coinciding borders of magnetic zone and there 
is complete nesting — flat parts of the Fermi surface coincide after translation by vector 
of antiferromagnetism Q = (±7r/a, ±7r/a). For T = electronic spectrum is unstable - 
energy gap is opened everywhere at the Fermi surface and the system becomes insulator due 
to formation of antiferromagnetic spin density wave (SDW) 9 These arguments lay behind the 



popular explanation of antiferromagnetism of HTSC - cuprates, let us quote e.g. Ref.[55 



The review of this and similar models can be found in 56 . In the case of the Fermi surface 



shown in Fig.[16] the appearance of antiferromagnetic long - range order, in accordance with 



general rules of the band theory ||57|| , leads to the appearance of discontinuities of isoenergetic 



surfaces (e.g. Fermi surface) at crossing points with borders of new (magnetic) Brillouin zone 



8 For example, tunneling experiments of Refs.|j2l[ 53] , in our opinion, practically rule out "superconducting" 
scenario of pseudogap formation. 

9 The same picture is valid for dielectrization due to formation of appropriate (period doubling) charge 
density wave (CDW). 
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FIG. 16: Fermi surface in the Brillouin zone and "hot spots" model. Dashed lines represent borders 
of magnetic Brillouin zone appearing after period doubling due to appearance of antiferromagnetic 
long - range order. "Hot spots" are defined as points of crossing of the Fermi surface with borders 
of magnetic zone. 

due to gap opening at points connected by vector Q. 

In the most part of underdoped region of the phase diagram of HTSC - cuprates there 
is no antiferromagnetic long - range order, however quite a number of experiments indi- 
cates the existence (below the T* line) of well developed fluctuations of antiferromagnetic 
short - range order. In the model of "nearly antiferromagnetic" Fermi liquid [[5^, |59 



an 



effective interaction of electrons with spin fluctuations is described by dynamic spin suscep- 
tibility x q (u;), which can be written in the following form determined from the fit to NMR 
experiments [j46 , 60|1 : 



K//(q, u) = g\^) - 1 + e2(q i Q)2 _._g_ (6) 

where g — is the coupling constant, £ - correlation length of spin fluctuations (short - range 
order), Q = (±7r/o, ±7r/a) - vector of antiferromagnetic ordering in insulating phase, u a f - 
characteristic frequency of spin fluctuations. 

Due to the fact that dynamic spin susceptibility Xqj\ u ) possesses peaks at wave vectors 
smeared around (±7r/a, ±7r/a), there appear "two types" of quasiparticles - "hot quasipar- 
ticles" with momenta from the vicinity of "hot spots" at the Fermi surface (Fig.|TBp and with 
energies satisfying the following inequality (vp - velocity at the Fermi surface): 

kk - £k+q| < v F /£, (7) 
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and "cold" ones, with momenta close to the parts of the Fermi surface around Brillouin 
zone diagonals \p x \ = \p y \ and not satisfying (^) |fT4"| . This terminology is due to the fact 
that quasiparticles from the vicinity of "hot spots" are strongly scattered with momentum 
transfer of the order of Q, because of strong interaction with spin fluctuations (B|), while for 
quasiparticles with momenta far from "hot spots" this interaction is much weaker. In the 
following this model will be called "hot spot" model 10 . 

Correlation length of short - range order fluctuations £, defined in (|6]), is very important 
for the following discussion. Note that in real HTSC - systems it is usually not very large 



and changes in the interval of 2a < £ < 8a [pi 



Characteristic frequency of spin fluctuations ou s f, depending on the compound and doping, 
is usually in the interval of 10 — lOOfT |61], |62[ , so that in most part of pseudogap region on 
the phase diagram we can satisfy an inequality 7rT ^> cu s f, which allows us to neglect spin 
dynamics and use quasistatic approximation 11 : 



-2 



Wq) = ^ 1 + p f q _ Q) , («) 

where W - is an effective parameter of dimensions of energy, which in the model of AFM 
fluctuations can be written as : 



W 2 = g 2 T £ x^uim) = g 2 <S 2 >/3 (9) 



mq 



where S« is the spin on the lattice site (ion of Cu in the CuO<i plane of HTSC - cuprates). 

In the following we shall consider both W and £ as phenomenological parameters. In 
particular, the value of W determines effective width of the pseudogap. The full micro- 
scopic theory of the pseudogap state will not be our aim, we shall only try to model the 
appropriate renormalization of electronic spectrum and its influence on physical properties, 
e.g. superconductivity. 

Significant simplification of calculations appears if instead of (^|) we use the model inter- 



10 Let us stress that the AFM nature of fluctuations is in fact irrelevant for further analysis and the model 
used in the following is, in fact, more general. The only thing important is the presence of strong scattering 
of electrons with momentum transfer from some narrow vicinity of the vector Q, which scatters electrons 
from one side of the Fermi surface to other. This type of scattering may be also due to charge (CDW) or 
structural fluctuations. 

11 In terms of Ref. fl4| this corresponds to the region of "weak" pseudogap. 
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action of the form (similar simplification was first used in Ref. ||ll|| ): 

wo = ^X-Q.n- + V-w (10) 

where W 2 = W 2 /4. In fact ( |10|) is qualitatively quite similar to (H) and differs from it only 
slightly quantitatively in most interesting region of |q — Q| < Analogous but a little bit 
different form of effective interaction was used in Ref . |I3J . In this way we achieve effectively 
one - dimensional situation. 

Scattering by AFM fluctuations in HTSC - cuprates is not always most intensive at 
vector Q = (it /a, vr/a) commensurate with the inverse lattice vectors, in general case vector 
Q may correspond to incommensurate scattering. Keeping in mind the observed topology 



of the Fermi surface with flat parts (patches) as shown in Fig.|12|, we may propose another 
model for scattering by fluctuations of AFM short - range order 12 , which we shall call the 
"hot patches" model [Q. Assume that the Fermi surface of our two - dimensional system 
of electrons has the form, shown in Fig.|I7]. The size of "hot patch" is determined by the 
angular parameter a. It is well known that the presence of flat parts on the Fermi surface 
leads to instabilities either of charge density wave (CDW) or spin density wave (SDW) type, 
with the establishment of the appropriate long - range order and opening of (dielectric) 
energy gap on these parts of the Fermi surface. Here we are interested in situation typical 
of fluctuation region of the appropriate phase transition, when there is still no long - range 
order present 13 . Similar model of the Fermi surface was introduced a number of years ago 



as a possible model for HTSC - cuprates in Refs. |]65|, p6[ p7|l , where rather detailed analysis 
of microscopic criteria for the appearance of antiferromagnetic and superconducting phases 
on the phase diagram was performed. Fluctuations of short - range order are assumed to 
be static and Gaussian 14 and the appropriate correlation function is defined as (similar to 



12 As in the following we always neglect, for simplicity, the spin structure of interaction, our analysis is, 
strictly speaking, applicable to the case of electron interaction with fluctuations of charge density wave 
(CDW) short - range order. However, this simplification is unimportant if we are interested only in 
qualitative results for pseudogap renormalization of electronic spectrum. Note that interaction with 
fluctuations of CDW type was proposed as an alternative mechanism of pseudogap anomalies of cuprates 
in a number of papers, see e.g.|Q. 

13 Assumption about the existence of flat parts on the Fermi surface is in fact not crucial for our model, but 

simplifies calculations considerably, which can always be done in more realistic "hot spot" model. 

14 This assumption, strictly speaking, is also applicable only for high enough temperatures. 
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FIG. 17: Fermi surface in the model of "hot patches", which are shown by thick lines (e.g. of the 



width of ~ £ ). The angle a determines the size of a "hot patch" [B3|, a = 7r/4 corresponds to 
the square Fermi surface. 



m 

5(q) = ^(g,-Q,) 2 + £- 2 (q y - Q y ) 2 + £ -2 (H) 

Here £ - is again the correlation length of fluctuations, while the scattering vector is taken 
as Q x = ±2kp, Q y = or Q y = ±2kp, Q x = 0. It is assumed that only electrons from flat 
( "hot" ) parts of the Fermi surface shown in Fig.[l7| interact with these fluctuations and this 
scattering is in fact one - dimensional. For the special case of a = 7r/4 we obtain square 
Fermi surface and the problem becomes strictly one - dimensional. For a < 7r/4 there 
are also "cold" parts present at the Fermi surface and the scattering there is either absent 
or weak enough. Effective interaction of electrons from "hot" patches with fluctuations 
will be descibed by (27r) 2 I4 y2 S'(q), where W is the parameter with dimension of energy, 
denning again the energy scale (width) of the pseudogap. 15 . The choice of scattering vector 
Q = (±2kp, 0) or Q = (0, ±2kp) assumes the picture of incommensurate fluctuations as 



15 It can be said that we introduce an effective interaction "constant" with fluctuations of the following form: 
W p = W[9( P x -p x )9{ P 'i+ Px ) + 9(p y ~p y )9( P Q y + Py )} (see Fig©. 
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the Fermi momentum pp = Ukp is in general not connected with the inverse lattice period. 
Commensurate case can also be studied within this model |k| . 

The basic idea of models under discussion is the presence of strong scattering on fluctu- 
ations of short - range order, which, according to (^), ( 110 ) or (|TTD, is effective in a limited 
region of inverse space with the width of the order of £ _1 around the "hot" spots or patches, 
and which leads to pseudogap - like renormalization of electronic spectrum in these regions 16 . 

Let us present a simple qualitative analysis of possible changes of single - electron spectral 
density (|5|). In the standard Fermi - liquid theory the single - electron Green's function 
for the metal can be written as: 

G{u, k) = f k , + G incoh (12) 

where ^ = — fi - is the quasiparticle energy with respect to the Fermi level (chemical 
potential) fj,, 7 k - its damping. The residue in the pole < Z k < 1, G inco h - non - singular 
contribution of many - particle excitations. Then the appropriate spectral density is: 



A(u>, k) = -Z k - - + ... (13) 



where dots denote more or less regular contribution due to Gi nco h, while the quasiparticle 
spectrum is reflected in a narrow (as damping 7k is small close to the Fermi level) Lorentzian 
peak. This situation is illustrated in Fig.|l8|(a). 

If the long - range order of SDW or CDW type appears the spectrum of elementary 
excitations acquires the (dielectric) gap (dependence on k stresses the possibility of gap 
opening on some part of the Fermi surface only) and single - electron Green's function has 



the form of Gorkov's function [55, 56 , where we can add some damping r k : 



G(u,k) = — — + — — (14) 



where the spectrum of elementary excitations: 



E* = Ja + W£ (15) 



Analogous situation is realized in liquid metals or semiconductors, where the information about the lost 
crystal structure is partially conserved in the so called structure factor with characteristic maximum in 
momentum space, being the crucial element of electronic theory of these systems according to Ziman and 
Edwards ]6§| . However, in three - dimensional isotropic liquids the angular averaging leads, in general, to 
rather strong suppression of pseudogap effects. 
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FIG. 18: Qualitative evolution of the spectral density, (a) - normal metal (Fermi - liquid), £k = 
- at the Fermi surface, (b) - two narrow peaks corresponding to Bogoliubov's quasiparticles in 
the system with dielectric gap (CDW or SDW long - range order). Smeared maxima - system 
without long - range order (pseudogap behavior), £k = - at the Fermi surface, (c) The same as 
in (b), but for £k > 0, i.e. above the Fermi surface. In this case a characteristic asymmetry of the 
spectral density appear. 

and we introduced Bogoliubov's coefficients: 

-* = 5 f 1 + t) (16) 

Then the spectral density: 



where now two peaks are present (both narrow if T k is small), corresponding to Bogoliubov's 
quasiparticles. 
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Let us now consider the situation when SDW or CDW long - range order is absent, 
but there is a strong scattering on fluctuations of short - range order with appropriate 
wave vectors (cf. (p|),([T0|) or (0)). Then it is not difficult to imagine that in the spectral 
density remains some "reminiscence" (or "anticipation") of dielectric energy gap opening 
(in the same region of momentum space) in the form of characteristic smeared two - peak 
structure, as shown in Fig.|TB| 17 . The width of both maxima, naturally, is determined by 
some parameter like vf/£, i.e. by inverse time of flight of an electron through the region 
of the typical size of £, where effectively persists "dielectric" ordering. It is seen that this 



qualitative picture is quite similar to ARPES data shown in Fig.[L5]. Below we shall see that 
more rigorous analysis leads just to the same conclusions. Further parts of the review will 
be dealing with justification of this qualitative picture and some of its consequences. 

Recurrence Relations for the Green's Functions. 

The model of the pseudogap state which will be discussed below is, in fact, a two - 
dimensional generalization of the model proposed many years ago for pseudogap formation 
in fluctuation region of Peierls (CDW) structural transition in one - dimensional systems 



ffO|) 71, 72, 73 1 . In Refs. fl7T|, [72| we obtained an exact solution of this model in the 



limit of very large correlation lengths of fluctuations of short - range order £ — ► oo, while 



in Ref. [73] we proposed the appropriate generalization for the case of finite £. The main 



advantage of these results, in comparison with Ref. [[70], was the complete account of all 
Feynman diagrams of perturbation series for electron interaction with fluctuations of short 
- range order. Appropriate two - dimensional generalization of this approach was made in 



Refs.JTJ], |I5J for "hot spots" model, while for "hot patches" model it can be done directly 
due to effective one - dimensionality of this model. 

Contribution of an arbitrary diagram for electron self - energy of single - particle Green's 
function in the N-th order of perturbation theory over interaction with fluctuations, given 
by ( |T0"D or (|ITD, can be approximately written in the following form [ I5 |: 

2JV-1 -i 

s (Af W) = w 2N n -. — r—. — (19) 



17 In the following the term "non - Fermi liquid behavior" will be used only in this, rather limited, sense. 

18 In Ref. this model was proposed for the explanation of pseudogap in liquid semiconductors 
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where for the "hot spots" model £j = £ p +q and Vj = \Vp +cl \ + [^p+qI for odd j and £j = £ p 
and Vj = \Vp\ + \v^\ for even j. In the "hot patches" model £j = (— 1) J £ P , = t>F. In (|I5| ) we 
have introduced an inverse correlation length k = - is the number of interaction lines, 

surrounding j-th (from the beginning) Green's function in a given diagram, e n = (2n + 1)ttT 
and, for definiteness, e n > 0. 

As a matter of fact, expression ([TJ]) constitutes rather successful Ansatz, which allows to 
calculate the contribution of an arbitrary diagram in any order. In this case any diagram 
with intersections of interaction lines is actually equal to some specific diagram of the same 
order without intersecting interaction lines. Thus, in fact we may consider only diagrams 
without intersections of interaction lines, taking those with intersections into account by 
additional combinatorial factors, associated with interaction lines. This approach was first 



used (for another problem) by Elyutin [74 and applied to one - dimensional model of the 



pseudogap state in Ref . |73[| . Actually, Ansatz (|T^) is not exact [[T^j, but it may be justified in 



two - dimensional case for certain topologies of the Fermi surface, such that electron velocity 
components in "hot spots" connected by vector Q have the same sign, so that fpfp+q > 



and v^Vp + Q > [15]. This condition is not satisfied in one - dimensional case and also for 
topologies of Fermi surfaces typical for HTSC - cuprates, as shown in Fig.|16|. However, in 
Ref. ]T3] it was shown that even in these cases the use of ( |T§D gives quite satisfactory results 
and reproduce the exactly known limits of £ —>■ oo and £ — > 19 . 

As a result we obtain the following recursion relation for one - electron Green's function 
(continuous chain approximation f73j): 

G-\e n £ p ) = G \e n £ p ) - £i(e n £ p ) (20) 

MenQ = W", V{k) (21) 

ISn ~ ffc + ikv k K - S fc+ i(e n 4 P ) 

Combinatorial factor: 

v(k) = k (22) 



19 Even in most unfavorable case of one - dimensional model Ansatz ( |l^ ) gives very good quantitative 
approximation, e.g. for the density of states, which is seen from the direct comparison f76fl of the results 
of |n| for incommensurate case with the results of an exact numerical simulation of Refs.j77], [7£|. In 
commensurate case the use of (|l|) in one - dimensional model does not reproduce only the so called 
Dyson singularity of the density of states at the center of pseudogap 0, f78| , which is apparently absent 
in two dimensions. 
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corresponds to the case of commensurate fluctuations with Q = (±n/a, ±n/a) fl73f . There 
is no problem to analyze also the case of incommensurate fluctuations when Q is not related 
to the inverse lattice period. In this case diagrams with interaction lines surrounding odd 
number of vertices are much smaller than those with interaction lines surrounding the even 
number of vertices. So we have to take into account these latter diagrams [ 71] , [72] , 75 . In 



these case the recursion relation Q2"T| ) is still valid, but diagram combinatorics change, so 
that combinatorial factor s v(k) are: |73| : 



. *±i for odd k 

v(k) = { (23) 

2 for even k 

In Ref . [[14]] the spin structure of interaction was taken into account in the framework of the 
model of "nearly antiferromagnetic" Fermi - liquid (spin - fermion model). This leads to 
more complicated diagram combinatorics for commensurate case with Q = (±7r/a, ±7r/a). 
Spin - conserving scattering leads formally to commensurate combinatorics, while spin - flip 
scattering is described by diagrams as for incommensurate case ( "charged" random field in 



terms of Ref.|T3|]). As a result, the recursion relation for the Green's function is still given 



by (pip, but the combinatorial factor v (k) becomes [I4J: 



, *±2 for odd A; 
v(k) = { I (24) 
for even k 
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Spectral Density and the Density of States. 

Recurrence relations ( p0|) and ( pT|) are quite convenient for numerical calculations. Rather 
detailed studies of spectral density, ARPES characteristics and density of states for different 
variants of "hot spots" model were performed in Refs. E , As a typical example in Fig.^9 



we present the results of |T5|] for the spectral density for incommensurate case. We can see 
that the spectral density close to the "hot spot" has the expected non - Fermi liquid like 
form, quasiparticles can not be defined. Far from the "hot spot" the spectral density reduces 
to a sharp peak, corresponding to well defined quasiparticles (Fermi - liquid). In Fig.^ from 



Ref. []14j we show the product of Fermi distribution and spectral density for different points 
on "renormalized" Fermi surface, defined by the equation — ReY>(uj = Ok) = 0, where the 
"bare" spectrum was taken as in (|3|) with t = — 0.25eV, t' = — 0.35t for hole concentration 
n h = 0.16, coupling constant in (|J) g = O.SeV^ and correlation length £ = 3a (commensurate 
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FIG. 19: Spectral density in "hot spots" model, incommensurate case, hm = 0.01 ]15|. (1) - at the 
"hot spot" PxCl/tt = 0.142, p y a/ir = 0.587. (2) - close to the "hot spot" for p x a/-K = 0.145, Pya/ir = 
0.843. (3) - far from the "hot spot" for p x a/n = p y /ir = 0.375. 

case of spin - fermion model). The qualitative agreement with ARPES data discussed above 
is clearly seen, as well as qualitatively different behavior close and far from the "hot spot" . 
In Fig.|21] from [14] we show calculated (for spin - fermion model with interaction static 
limit) positions of the maximum of A(uh) for two different concentrations of holes and 
appropriate experimental ARPES - data of Ref.|79[[ for Bi 2 Sr2Cai- x Dy x Cu20s+s- These 



maxima positions in (u, k) plane, obtained from ARPES, in fact determine (for Fermi - 
liquid like system) the experimental dispersion curves of quasiparticles (Cf. Fig.|Tj|(a)). For 
overdoped system the values of nh = 0.22 and £ = a were assumed. The results show rather 
well defined branches of the spectrum both in the direction of diagonal of the Brillouin zone, 
as well as in the (0,0) — (7r, 0) direction. For underdoped system the values of rih = 0.16 
and £ = 3a were used. In this case, in diagonal direction we can again see the crossing of 
the spectrum with the Fermi level, while close to the "hot spots" in the vicinity of (it, 0) 
the smeared maximum of the spectral density remains approximately 200meV lower than 
the Fermi level (pseudogap behavior). In general, the agreement of theoretical model and 
experiment is rather encouraging. 
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FIG. 20: Product of spectral density and Fermi distribution function at different points of the 
Fermi surface, shown at the insert. Spin - fermion model, correlation length £ = 3a |L4| l. 
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FIG. 21: Position of the maximum of spectral density for two values of correlation length £ and 
hole concentration, calculated for spin - fermion model |l4|] in comparison with ARPES data of 
Ref.]79[| for Bi2Sr2Ca\- x Dy x Cu20% + $ with x = 1 (triangles) and x = 0.175 (diamonds). 
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Consider now the single - electron density of states: 



N(E) = £ A(E, P ) = -- E ImG R (Ep) (25) 



p * p 



defined by the integral of spectral density A(Ep) over Brillouin zone. Detailed calculations 
of the density of states in the "hot spots" model were performed in Ref.||15||. For typical 
value of t'/t = —OA rather smooth suppression of the density of states (pseudogap) was 
obtained. This drop in the density of states is relatively weakly dependent on the value 
of correlation length £. At the same time, for t'/t = —0.6 (non characteristic for HTSC - 
cuprates) the "hot spots" on the Fermi surface still exist, but pseudogap in the density of 
states is practically invisible. Only the strong smearing of Van - Hove singularity, existing 
for the system without fluctuations, is clearly observed. 

In "hot patches" model the use of (p0|) and (pID leads to the spectral density on these 
patches qualitatively similar to that shown in Fig.[l9| |8l|, jB2|j . On "cold" parts of the Fermi 
surface the spectral density reduces to ^-function, typical for Fermi - liquid (gas) (Cf. Fig.|l8 
(a)). The density of states can be written as: 

oo 

N(E) = --iVo(O) r ^ / d{ p ImG R (EQ = —N W {E) + (1 - — )N (0) (26) 
n Jo 2n J 7r 7r 

— oo 

where N (0) - is the density of states of free electrons at the Fermi level, while N\y(E) - is 
the density of states of one - dimensional model (square Fermi surface) , obtained previously 



in Refs. 71, 72, 73 



Let us quote some more detailed results for (artificial enough) limit of very large corre- 
lation lengths £ — > oo ("toy" - model). In this case the whole perturbation series for the 



electron scattered by fluctuations is easily summed [71, 72] and an exact analytic solution 



for the Green' function is obtained in the following form [33]: 



O^P) = JdDnD) ^ (27) 



where D((j)) is defined for < < | as: 



, D ,0<<f><a, §-«<</><§ 
D{<j>) = { ~ ~ 2 ~ ~ 2 (28) 



,a < d) < f -a 



For other values of the value of D(4>) is determined by obvious symmetry analogously to 

34 



The amplitude of dielectric gap D is random and distributed according to Rayleigh distri- 
bution ]7vJ (its phase is also random and distributed homogeneously on the interval (0, 2%)): 

2D ( D 2 \ 
T{D) = — exp {--) (29) 

Thus, on the "hot patches" the Green's function has the form of "normal" Gorkov's function, 
averaged over spatially homogeneous fluctuations of dielectric gap D, distributed according 



to [29]. "Anomalous" Gorkov's function is, however, zero due to random phases of dielectric 
gap D, which corresponds to the absence of the long - range order 20 . 

For finite correlation lengths £ the amplitude of our one - dimensional random "periodic" 
field is approximately constant on the lengths of the order of £, while its values in different 
regions with sizes of the order of ~ £ are random. The crossover from one region with 
approximately periodic field to the other also takes place at lengths of the order of £. The 
electron is effectively scattered only during this transition from one region to another, which 
takes a time of the order of ~ £/vp, leading to effective damping of the order of ~ vp/£. 
Interesting data on numerical modeling of our random field, for one - dimensional case, can 



be found in Ref. |75| . 

Outside "hot patches" (the second inequality in fl2lf)) the Green's function (E7I) just 
coincides with that of free electrons. 

Density of states, corresponding to (^7j), has the form (|26|) , where 



N W {E) 
iVo(O) 



E 



W 



o y w 2 ^ 



where Erfi(x) - is the probability integral of imaginary argument. In Fig .^2] we show 
density of states of our model for different values of parameter a, i.e. for "hot patches" of 
different sizes. It is seen that the pseudogap in the density of states is rather fast smeared 
(or filled in) as "hot patches" become smaller. In some sense the effect of diminishing a is 
analogous to the effect of diminishing correlation length £ J73 |, so that the approximation 



of £ — > oo, is probably not so very strong limitation of this model. Finite values of £ can 
be easily accounted using (EH) and (^T|), leading to additional smearing of the pseudogap 



20 Note that averages of pairs of anomalous functions are different from zero and contribute to the appro- 
priate exact expression for two - particle Green's function jnj [T2| 
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FIG. 22: Density of states for "hot patches" of different sizes |63|1 : (1) - ol = vr/4; (2) - a = 7r/6; 
(3) - a = vr/8; (4) - a = vr/12; (5) - a = vr/24. 

for smaller £. Note the general qualitative agreement of the form of the pseudogap in this 
model with those observed in tunneling experiments, e.g. shown in Figs.^J^. 

Two-Particle Green' Function and Optical Conductivity. 

Remarkably, for these models it is possible to calculate explicitly also two - particle 
Green's functions for the electron in the random field of short - range order fluctuations 



81] (Cf. also Ref. [0), taking into account all Feynman diagrams of perturbation 
series. As calculations for "hot spots" model require rather large numerical work due to the 
use of "realistic" spectrum of current carriers (|3]), we shall limit ourselves to the simplified 



analysis within "hot patches" model |83 



On "cold" parts of the Fermi surface we shall assume the existence of some small static 
scattering of an arbitrary nature, describing it by phenomenological scattering frequency 
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7, assuming, in most cases, that 7 <C W, so that we may neglect this scattering on "hot 
patches" . Accordingly, on "cold" parts the electronic spectrum will be described by by usual 
expressions for the Green's functions for the system with weak (static) scattering. 

Consider first the limit of £ —>■ 00 , when the single - electron Green's function has the 
form ( P7D and two - particle Green's function can also be found exactly using methods of 
Refs. J7I|,|723. 

Conductivity in this model always consists of additive contributions from "hot" and 
"cold" parts of the Fermi surface, analogously to (26). In particular, for its real part in the 
limit of £ —> 00 we have: 



4a 4a 
Rea(uj) = — Reaw(u) + (1 )Rea£>{uj) 

IX 7T 



where, with the account of the results of Refs. [[Tl|, [72]] : 

2 ^ u 2 /AW 2 

Rea w (uj) = / d(exp(-() 

4 uj z J 


where uj p - is plasma frequency and 

Rear, 



c 



uj 2 /4W 2 -( 



7 



(31) 



(32) 



(33) 



Aix uj 2 + 7 2 

- the usual Drude - like contribution of "cold" parts. 

Even in this simplest approximation the dependence of Rea on uj is quite similar to 

]. It is characterized by narrow "Drude" peak at 



observed in the experiments |29], |3(J, [84], 
small frequencies and smooth pseudogap absorption maximum at uj ~ 2W. As scattering 
frequency 7 on "cold" parts grows, characteristic "Drude" peak at small frequencies is 
suppressed. 

More realistic case of finite correlation lengths of fluctuations of "antiferromagnetic" 
short range order £ in ([11]) can be analyzed using (go)), (|21|). Vertex part J RA (e,£ p ;e + 
a;,£ p+q ), defining density - density response function (two - particle Green' function) on 
"hot patches" an be determined from the following recursion procedure (details can be found 
in Refs. 1 



JTJ and also in |I4"|), which takes into account all diagrams of perturbation series 
for interaction with fluctuations: 



x 1 + 



Jk-l( 6 l 6 + W i £p+q) - 

e + W 2 v(k)G£(e, £ p )Gf (e + uj, £ p ; e + u, £ p+q ) > 

Hvp nk 

uj - (-l) k v F q + v(k + l)W 2 [G£ +1 (e, £ p ) - Gf +1 (e + uj, £ p+q )] 
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FIG. 23: Real part of conductivity as a function of frequency of external field, for different values of 
7 and fixed correlation length vpK = 0.5W. Conductivity is in units of LOp/AnW. (l)—j/W = 0.2; 
(2)-~f/W = 0.5; (3)-7/W = 1.0. The size of "hot patches" a = vr/6. §||. At the insert - calculated 
dependence of generalized scattering frequency 1/r (in units of W), as a function of frequency for 
different values of correlation length and fixed 7 = 0.2W. (1)-vfk = 0.1; (2)-vfk = 0.5; (3)- 



= 1.0 mi 



where e - is electron charge, R(A) denotes retarded (advanced) Green's function, and the 
vertex part of interest to us is denned as k = term of this recursion procedure. Contribution 
of "hot patches" to conductivity Rea w {u) in ([H]) is calculated as in Refs. pC| |81| , while 
Reaoioj) is again determined by (|33|). 



Typical results of such calculations are shown in Fig.[23| for the case of incommensurate 
fluctuations. The use of combinatorics of spin - fermion model leads to rather insignificant 
quantitative changes. The general qualitative picture remains the same also for commensu- 
rate case ||831 . Conductivity is always characterized by narrow enough "Drude" peak in the 
region of small frequencies u < 7, originating from "cold" parts of the Fermi surface and a 
smooth maximum in the region of ui ~ 2W, corresponding to the absorption through the 
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pseudogap opened at the "hot patches'.' "Drude" peak is rapidly suppressed as 7 grows. 
Dependence on correlation length of fluctuations £ = n^ 1 in most interesting parameter 
region is weak enough. This qualitative picture is very similar to experimental dependencies 
observed for many HTSC - systems P9"| , Q |8"5| , typical data were shown above in Fig.|S|. 



In this model we can also calculate the imaginary part of conductivity and the parameters 
of the generalized Drude model ([TJ), (0) Calculated 1/t(u) is shown at the insert in 



Fig.|23| and demonstrates pseudogap behavior in the frequency region of lj < W, quite anal- 
ogous to that shown in Fig.|l(] 21 . Apparently it is not difficult to fit experimental data to 
our dependences using the known values of u p ~ 1.5 — 2.5eV and 2W ~ O.leV, as well as 
"experimental " values of 7, determined from the width of "Drude" peak, and varying "free" 
parameters a (the size of "hot patches") and £ (the value of this can also be taken from 



other experiments fH 



SUPERCONDUCTIVITY IN THE PSEUDOGAP STATE. 
Gorkov's Equations. 

At present there is no agreement on microscopic mechanism of Cooper pairing in HTSC 
- cuprates. The only well established fact is the anisotropic nature of pairing with d-wave 



symmetry f45| , f47|, though even here some alternative points of view exist. Apparently, 
the majority of researchers are inclined to believe in some kind pairing mechanism due 
to an exchange by spin (antiferromagnetic) excitations. A recent review of these theories 
can be found in Ref. ||46|| . Typical example of such a theory is the model of "nearly 
antiferromagnetic" Fermi - liquid (spin - fermion model), intensively developed by Pines and 



his collaborators [58 , [59]. This model assumes that electrons interact with spin fluctuations 



and the form of this interaction (|6T) is restored, as mentioned above, from the fit to NMR 



experiments. [50|. This approach continues to develop rapidly, we can only mention the 



very recent papers |S7L BR B3] . There exist rather convincing experimental facts supporting 



this mechanism of pairing H, pOj. At the same time interaction of the form of (o) can be 



also responsible, as we shown above, for the formation of "dielectric" pseudogap already at 



21 For u> > 214^ the linear growth of 1/t(lo) with frequency is not reproduced, but it can be easily obtained 
if we use a phenomenological dependence 7(0;) = 70 + aw in the spirit of "marginal" Fermi - liquid. 
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high enough temperatures. Unfortunately, up to now there are no papers analyzing both 
the pseudogap formation and superconductivity within the single approach of spin - fermion 
model. 

On the other hand rephrasing the quotation from Landau on Coulomb interaction, cited 
in |91|| , we can only say that electron - phonon interaction was also "not canceled by any- 
body" , as a matter of fact as a possible mechanism of pairing in cuprates also. The review of 
appropriate calculations with rather impressive results can be found in |92| and in earlier 



Ref . |93| . Probably the main difficulty for the models based on electron - phonon interaction 
in cuprates remains the established <i-wave symmetry of Cooper pairs. However, there exist 
interesting attempts to explain the anisotropic pairing within the electron - phonon mecha- 
nism [91], as well as serious doubts in the effectiveness of spin - fluctuation mechanism 

The author of the present review is not inclined to join either of these positions, being 



convinced that these (or similar [94]) more or less "traditional" mechanisms are apparently 
decisive for the microscopies of HTSC. Anyhow, in the following we shall not even discuss 
some other more "radical" approaches, e.g. those developed in different variants and for 
many years models, based on the concept of Luttinger liquid ||9~5| , mainly because of the 
absence there of any "stable" results and conclusions. 

Significant renormalization of electronic spectrum due to the pseudogap formation in- 
evitably leads to a strong change of properties of the system in superconducting state. In 
the scenario of pseudogap formation due to fluctuations of "dielectric" short - range order 
(AFM, SDW or CDW) we can speak precisely about the influence of the pseudogap upon 
superconductivity. In this formulation the analysis can be done even without assuming 
specific pairing mechanism, as it is done e.g. in the problem of the influence of structural 
disordering, impurities etc. From the physical point of view the least justified is the used 
above assumption of static nature of short - range order fluctuations, as dynamics of these 
fluctuations (e.g. in the model based on (|6])) is apparently decisive at low enough temper- 
atures (in superconducting phase), being responsible also for the pairing mechanism itself. 
We deliberately use these simplification as the full self - consistent solution of the problem 
within the model with dynamics seems to be out of reach, at least at present. Besides, in 
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the following we shall only use the most simplified "hot patches" model 22 . For generality 
we analyze both s-wave and d-wave pairing. 

So, in the spirit of the previous discussion we shall not concentrate on any concrete 
microscopic mechanism of pairing and just assume the simplest possible BCS - like model, 
where the pairing is due to some kind attractive "potential" of the following form: 

V sc (p, p') = V(<j>, 0') = -Ve(0)e(0'), (35) 

where - is polar angle defining the direction of electron momentum p in highly - conducting 
plane, while for e(0) we shall use the simple model dependence: 

e(0) = ( 1 ( Pairing) . (36) 

I -\/2cos(20) ( (i-wave pairing) 

Attractive interaction constant V is assumed, as usual, to be nonzero in some energy interval 
of the width 2u c around the Fermi level (u c - is some characteristic frequency of excitations, 
responsible for attraction of electrons). The model interaction of the type of ([35|) was 
successfully used e.g. in Refs.|)6], |97j for the analysis of impurities influence on anisotropic 
Cooper pairing. 

Superconducting energy gap in this case takes the following form: 

A(p)=A(0) = Ae(0). (37) 

In the following, just not to make expressions too long, we shall use the notation of A instead 
of A(0), writing the angular dependence explicitly only when it will be necessary. 

In superconducting state the perturbation theory for electron interaction with AFM fluc- 



tuations (|TTD should be built on "free" normal and anomalous Green's functions of a super- 
conductor: 

G 00 (e n p) = - 2 , « _5|A|2 ; F oo( £ nP) = 2 f2 |A|2 (38) 

n ' ' I I n ' ^p ' II 



In this case we can formulate the direct analog of approximation (|T9| ) also in superconducting 
phase |98 1 . The contribution of an arbitrary diagram of the 7V-th order in interaction (flO|) to 
either normal or anomalous Green' function consists of a product of A^+l "free" normal Gok 



22 This simplification is not very important. The same analysis can be performed also in more realistic "hot 
spots" model, but this leads to much more complicated (numerically) calculations. 
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or anomalous F^ k . Green' functions with renormalized (in a certain way) frequencies and 
gaps. Here kj - is the number of interaction lines surrounding the j-th (from the beginning 
of the diagram) electronic line. As in normal phase the contribution of an arbitrary diagram 
is defined by some set of integers kj, and each diagram with intersecting interaction lines 
becomes equal to some diagram without intersections. For this last reason we can analyze 
only diagrams with no intersections of interaction lines, taking all the other diagrams into 
account by the same combinatorial factors v(k), attributed to interaction lines as in the 



normal phase. As a result we obtain an analog of Gorkov's equations ||69fl . Accordingly 
there appears now a system of two recursion equations for normal and anomalous Green's 
functions: 

Gk = G k + GokGGk — G ok FF k — F ok G*F k — F ok F + G k 

F k = F ok + F okGG k - F ok FF k + G* Qk G*F£ + G* ok F + G k (39) 

where 

G = W 2 v(k + l)G*+i; F + = W 2 v{k + 1)F+ +1 (40) 

G ok (E n p) = _ z£ n + (- 1 )% F +( £nP ) = A * _ (41) 

and we introduced the abovementioned renormalized frequency e and gap A: 

~ kvpK 
rj k e n ; A = r] k A; r} k = l+ - (42) 

lei + |A| 2 



similar to those appearing in the analysis of a superconductor with impurities |69| . 

Both normal and anomalous Green's functions of a superconductor with pseudogap are 
defined from (|3^) at k — and represent the complete sum of all diagrams of perturbation 
series over interaction of superconducting electron with fluctuations of AFM short - range 
order. 

In fact we just consider Gorkov's Green's functions averaged over an ensemble of random 
(Gaussian) fluctuations of short - range order, in many respect analogously to the similar 
analysis of the case of impure superconductors |69j . Here we assume here the self - averaging 
property of superconducting order parameter (energy gap) A, which allows to average it 
independently of electronic Green's functions in diagrammatic series. Usual arguments for 



the possibility of such independent averaging are the following |M 100, 101]: the value 



of A changes on characteristic length scale of the order of coherence length £ ~ ^f/A of 
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BCS theory, while Green's functions oscillate on much shorter length scale of the order of 
interatomic distance. Naturally, this last assumption becomes wrong if a new characteristic 
length £ — > oo appears in electronic subsystem. At the same time, until this AFM correlation 
length £ <C £o (i- e - when AFM fluctuations correlate on lengths smaller than characteristic 
size of Cooper pairs), the assumption of self - averaging nature of A is apparently still 
conserved, being broken only in the region of £ > £o- As a result we can use the standard 
approach of the theory of disordered superconductors (mean - field approximation in terms 
of Ref.[ 102|| ). Possible role of non - self - averaging A [ 102|| will be analyzed later. Note that 



in real HTSC - cuprates apparently always we have £ ~ £o; so that they belong to the region 
of parameters most difficult for theory. 

Critical Temperature and the Temperature Dependence of the Gap. 

Energy gap of a superconductor is defined by the equation: 

A(p) = -T £ £ V sc (p, p')F(e nP ') (43) 

p' £n 

On flat parts of the Fermi surface the anomalous Green's function is determined by recur- 
sion procedure (|3^). On the rest ("cold") part of the Fermi surface scattering by AFM 
fluctuations is absent (in our model) and the anomalous Green's function has the form as 
in Eq . (|38|) . The results of numerical calculations of the temperature dependence of the gap 
for different values of correlation length of short - range order fluctuations £ can be found 
in Ref . ||| . These dependences are more or less of the usual form. 

Equation for the critical temperature of superconducting transition T c follows immedi- 
ately from (|43"D as A — > 0. Calculated in Ref.|98[| dependences of T c on the width of the 
pseudogap W and correlation length (parameter k = £ _1 ) are shown in Fig.|24], where T c q - is 



transition temperature in the absence of the pseudogap. In this calculations we have, rather 
arbitrarily, taken the value of a = 7r/6, which is pretty close to the experimental data of 
Ref. . The general qualitative conclusion is that pseudogap suppresses superconductivity 



due to partial "dielectrization" of electronic spectrum on "hot patches" of the Fermi sur- 
face. This suppression is maximal for k = (infinite correlation length of AFM fluctuations) 
f63| , |102|1 and diminishes with diminishing values of correlation length, which qualitatively 
corresponds to the observed phase diagram of HTSC - systems. As noted above, parameters 
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FIG. 24: Dependence of superconducting critical temperature on the width of the pseudogap W 
and correlation length of AFM fluctuations (parameter k = ^ = 0.1 (1); = 1.0 (2); 

Zffi = 10.0 (3). Dashed line — n = Jfr|. At the insert: dependence of T c on k for ^ = 5 @ . 



of our model are to be considered as phenomenological. E.g. the effective width of the 
pseudogap 2V4 7 may be, apparently, identified as E g , the experimental values of which and 
dependence on doping are shown in Fig.§ for YBCO. Experimental data on the values of 
correlation length £ and its temperature dependence are rather incomplete. Indirectly this 
information may be extracted from NMR experiments [|14"|, |60[1. Direct neutron scattering 



data are not conclusive enough. As an example we shall mention the paper [|103|| , where the 
authors present some compilation of data on effective width of neutron scattering peak at 
the vector (ir/a,7r/a) in YBCO with different oxygen content. The inverse width of this 



peak can be naturally identified with the value of £. The authors of Ref. ||103|| observed an 
interesting correlation of superconducting critical temperature and thus defined values of 
which is shown in Fig.|25]. Obviously, this dependence is in direct qualitative corre- 
spondence with dependence of T c on k = shown in Fig.^j. Quantitative fit of these 
data to calculated dependences, following from our model, is apparently possible, but we 
must be also aware of rather incomplete information on doping dependence of parameter 



W. In particular, it is not clear, whether there is any deep physical meaning of E g 



W 



going to zero at some "critical" doping as shown in fig.|| Pseudogap influence on electronic 
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FIG. 25: Temperature of superconducting transition T c as a function of the width of neutron 
scattering peak at q = (it /a, vr/a), which is identified with the inverse of correlation - length of 
short - range order fluctuations |103|]. 



properties may disappear also due to correlation length £ becoming just smaller, which leads 
not to the "closing", but to the "filling" of the pseudogap in the density of states. However, 
if we accept the data of Fig.||] on the width of the pseudogap |J, then our parameter 2W 
drops from the values of the order of 7007^ at hole concentration p = 0.05 to the values 
of the order of T c ~ 100i^ near optimal doping p = 0.17, going to zero at p = 0.19. If 
we assume the validity of microscopic expression (|9]) for W, following from the theory of 
"nearly antiferromagnetic" Fermi - liquid, this concentration dependence may be attributed 
to the appropriate dependence of local spin < Sf > on Cu - ion. The author is unaware of 
any direct data of this sort, but we may note Ref.| 104 |, where the disappearance of effective 
interaction with antiferromagnetic fluctuations at p — 0.19 was demonstrated from some fit 
to experimental data on NMR relaxation. 

Let us stress once again, that the above theoretical results are valid only assuming self - 
averaging nature of superconducting order parameter (gap) as averaged over AFM fluctua- 
tions (random mean - field approximation ||102|| ), which is apparently true for not very large 
values of correlation length £ < £o ; where £o ~ is the coherence length of a superconductor 
(the size of Cooper pairs at T = 0). For £ 3> £ , as we shall see below, important effects due 
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the absence of self - averaging appear, leading to characteristic "tails" in the temperature 
dependence of the average gap in the region of T c < T < T co [ |102| . 



Cooper Instability. 

It is well known that the critical temperature may be determined also in another way, 
that is from the equation for Cooper instability of the normal phase: 

l-y X (0,0) = (44) 

where x(0, 0) is the generalized Cooper susceptibility: 

X (q0; T) = -T ]T ]T G(e n p + q)G(-£ n p)e 2 (0)r(e n , -e n , q) (45) 

Here we have to calculate the "triangular" vertex part T(e n ,— e n , q), describing electron 
interaction with AFM fluctuations. For one - dimensional analogue of our problem (and for 
real frequencies, T = 0) appropriate recursion procedure was formulated in Refs.[[80]. For 
our two - dimensional model this procedure was used in calculations of optical conductivity, 
as discussed above p3(] . This procedure is easily generalized for the case of Matsubara 
frequencies P5| . Below, for defmiteness, we assume e n > 0. Then, similarly to (|H4]), we 
have: 

r fc _i(e n , -e n , q) = 1 + 

+W 2 v(k)G k G Jl + — - k ^ W„, s n , q) 

( 2ie n - {-l) k v F q - W 2 v{k + l)(G*+i - G k +i) J 

r(e n , -e n , q) = r (£ n , -e n , q) (46) 

where G^ = Gk(£ n P + q) and G^ = Gk(—£ n , p) are calculated according to (|2l"D . 

To calculate T c we need the vertex at q = 0. Then Gt = G* k and vertices Y k become real, 
which leads to considerable simplification of the procedure (|46|) . 

The following exact relation holds, similar to Ward identity f98| : 

G(e n p)G(-e n p)T(e n , -e B , 0) = J mG ^ (47) 

£-n 

Numerical calculations fully confirm this relation, demonstrating the self - consistency of our 
recurrence relations for single - particle Green's function and vertex part. The relation ( f47|) 
leads to the equivalence of T c equation, derived from Cooper instability, and T c equation, 
obtained from the linearized gap equation, though the recurrence procedures used to derive 
both equations (accounting for AFM fluctuations) seem to be quite different. 
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Ginzburg-Landau Equations and Basic Properties of a Superconductor with the 
Pseudogap close to T c . 



Ginzburg - Landau expansion for the difference of free - energy density of superconducting 
and normal state can be written in standard form: 

F s - F n = A|A q | 2 + g 2 C|A q | 2 + || A q | 4 , (48) 

where A q - is the amplitude of Fourier - component of superconducting order parameter: 

A(0,q) = A q e(0). (49) 

This expansion (|48"D is determined by loop - expansion diagrams for the free - energy of an 
electron in the field of random fluctuations of order - parameter with small wave - vector q 

H- 

Let us write the coefficients of Ginzburg - Landau expansion in the following form: 

A = A K A ; C = C K C ; B = B K B , (50) 

where Aq, Cq and Bq are standard expressions for these coefficients for the case of isotropic 
s-wave pairing: 

^o = iV (0)^; C = N (0) 7 ^± Bo = JVb (0)^, (51) 

The all changes of these coefficients due to the appearance of the pseudogap are contained 
in dimensionless coefficients Ka, Kc and Kb- In the absence of the pseudogap all these 
coefficients are equal to 1, except Kb = 3/2 in case of <i-wave pairing. 

Consider again the generalized Cooper susceptibility (|45|). Then the coefficients Ka and 
Kc can be written as PS : 



X(00;T)- X (00;r c ) 
a = 2 ( 52 J 

x(qO;T c )-x(00;T c ) 
K c = km ^ (53) 



Then all calculations can be made using the recurrence procedure 

To find the coefficient B, in general case, we need more complicated procedure. Significant 
simplifications appear if we limit ourselves, for terms of the order of |A q | 4 in free energy, 
by the usual case of q — 0. Then the coefficient B may be determined directly from the 
anomalous Green's function F, i.e. from the recurrence procedure ( |39D ||98|| . 
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In the limit of £ — > oo all coefficients of Ginzburg - Landau expansion may be obtained 
analytically p5fl , using the exact solution for the Green's functions discussed above. 



Ginzburg - Landau expansion defines two characteristic lengths of superconducting state: 
coherence length and penetration depth (of magnetic field). Coherence length at given 
temperature £(T) determines the length scale of inhomogeneities of the order parameter A: 



e(T) = ~. (54) 



In the absence of the pseudogap: 



e BC s(T) = (55) 

Zbcs(T) « 0.74 e ° , (56) 
! l - T/T c 



where £o = 0.18vf/T c . In our model: 

^ - K ° (57) 



Appropriate dependences of i 2 {T) /^ 2 BCS {T) on the width of the pseudogap W and correlation 
length of fluctuations (parameter k) , for the case of <i-wave pairing were calculated in Ref. BH 
and found to be smooth, with relatively insignificant changes of the ratio ( |57D itself. 
For the penetration depth of a superconductor without the pseudogap we have: 

X BCS (T) = i= , A ° , (58) 
yfty/l- T/T c 

where Aq = 4 ™^ 2 defines penetration depth at T = 0. In general case: 

A2 < T) - "SOT < 59 » 



Then for our model we have: 



A < T > -I K ° V /2 . (60) 



\ BCS (T) \K a K c . 

Dependences of this ratio on the pseudogap width and correlation length were also calculated 
for the case of (i-wave pairing in Ref.[]98| and found to be relatively small. 



as: 



Close to T c the upper critical field H c2 is determined via Ginzburg - Landau coefficients 

0o _ 0o A 
Hc2 ~ 2^HT) ~ "2^C' (61) 
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where 0o = cir/e - is the magnetic flux quantum. Then the slope of the upper critical field 
near T c is given by: 



dH, 



<-2 



dT 



24tt0o „ K A 

- T c~r^- ( 62 ) 



7C(3)4 C K C 

Graphic dependences of this slope ^f- ^ , at temperature T c0 on the effective width of the 

pseudogap W and correlation length parameter k can also be found in Ref.|P5|. The slope 

of the upper critical field for the case of large enough correlation lengths drops fast with the 

growth of the pseudogap width. However, for short enough correlation lengths some growth 

of this slope can also be observed at small pseudogap widths. At fixed pseudogap width the 

slope of H c2 grows significantly with the drop of correlation length of fluctuations. 

Finally, consider the specific heat discontinuity (jump) at the transition point: 

C a C n T r ( A \ . . 

' (63) 



to B \T — T, 



- c 



where C s , C n - are the specific heats in superconducting and normal states respectively, 
Q, - is the sample volume. At temperature T = T c0 and in the absence of the pseudogap 
(W = 0): 



C s -C n \ ^^SttX 



(65) 



n ^ "<°>W (64) 

Then the relative specific heat discontinuity in our model can be written as: 

(C s — C n )T c T c K\ 
[C s — C n )T c0 T c0 Kb 

Appropriate dependences on effective width of the pseudogap W and correlation length 
parameter k for the case of <i-wave pairing are shown in Fig.|26]. It is seen that the specific heat 
discontinuity is rapidly suppressed with the growth of the pseudogap width and, inversely, 
it grows with for shorter correlation lengths of AFM fluctuations. 

For superconductors with s-wave pairing all the dependences are quite similar to that 
found above, the only difference is the larger scale of W leading to significant changes of 
all characteristics of superconducting state, corresponding to larger stability of isotropic 
superconductors to partial "dielectrization" of electronic spectrum due to to pseudogap 



formation on "hot patches" of the Fermi surface |63| , |102| . 

All these results are in complete qualitative agreement with experimental data on specific 
heat discontinuity discussed above @, 0, 0. As can be seen from Fig.||, specific heat dis- 
continuity is rapidly suppressed as system moves to underdoped region, where the pseudogap 
width grows as is seen from the data shown in Fig. 
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FIG. 26: Dependence of specific heat discontinuity on pseudogap width W and correlation length of 



AFM fluctuations (parameter k = f 1 ): Hffi = 0.1 (1); ^ = 1.0 (2); ^ = 10.0 (3). Dashed 



w 

line — k = p3] . At the insert: dependence of specific heat discontinuity on k for ^ 



5. 



Effects due to non-Self- Averaging Nature of the Order-Parameter. 



The above analysis of superconductivity assumed the self - averaging property of super- 
conducting order parameter A. This assumption is apparently justified in case of correlation 
length of fluctuations of (AFM) short - range order £ being small in comparison with char- 
acteristic size of Cooper pairs £o (coherence length of BCS theory). The opposite limit 
of £ 3> £o can be analyzed within the exactly solvable model of the pseudogap state with 
£ — > oo, which was described above for the case of "hot patches" model in (p7[) , (|29|) ||102|| . 

Consider first superconductivity in a system with fixed dielectric gap D on "hot patches" 
of the Fermi surface. The problem of superconductivity in a system with partial dielec- 
trization of electron energy spectrum on certain parts of the Fermi surface was considered 
previously in a number of papers (cf. e.g. Refs. ||105| , |106|| ), and in the form most useful for 
our case in a paper by Bilbro and McMillan |[L07|| . Below we shall use some of the results 
of this paper. Again we assume the simplest BCS - like model of pairing interaction (p5|). 



For the fixed value of dielectric gap D on "hot patches" of the Fermi surface supercon- 
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ducting gap equation, defining A for the case of s-wave pairing, takes the following form: 
A / d£{a j= _ 2T + (l-a) — 2T } (66) 



o 



where A = ViVo(O) - is dimensionless coupling constant of pairing interaction, a — — . The 
first term in the r.h.s. of Eq.(|66|) corresponds to the contribution of "hot" (dielectrized or 
"closed") parts of the Fermi surface, where electronic spectrum is given by ||107|| : E p = 



Q + D 2 + A 2 , while the second term gives the contribution of "cold" (metallic or "open") 



parts, where the spectrum has the usual BCS form: E p = y\£ 2 + A 2 . Eq. (p6|) determines 
superconducting gap A(D) at fixed value of dielectric gap D, which is different from zero 
on "hot patches". 

In case of <i-wave pairing the analogous equation has the form: 

1 = A- / dd #e 2 (0) : 2T + / d<Pe 2 (<P)- 2T 



""■>" I 7 11 J£ 2 + D 2 + A 2 (D)e 2 (0) •><» v /^ 2 + A 2 (D)e 2 (0) 

(67) 

From these equation it can be seen that A(D) diminishes with the growth of D, while A(0) 
coincides with the usual gap Ao in the absence of dielectrization on flat parts of the Fermi 
surface and appearing at temperature T = T c0 , defined by the equation: 

1 = A / d£— ^ (68) 
Jo 4 

both for s-wave and <i-wave pairing. 

For D — > oo the first terms of Eqs.(|66|),(|6"T|) tend to zero, so that the appropriate equations 
for Aqo = A(D — > oo) are: 



A / df(l — a) — 2T (s-wave pairing) (69) 
Jo ' & J- A2 



4 r«c r /4 th Ve±^M 
A- / d£ / ddx?{6)—= 2T (rf-wave pairing) (70) 

Wo L J^ + A^e 2 ^) 



Equation ( |6TSD coincides with gap equation for the case of D = with "renormalized" 
coupling constant A = A(l — a), so that for the case of s-wave pairing: 

A 00 = A (A = A(l-d)) (71) 
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and non zero superconducting gap in the limit of D — > oo appears at T < T coo 

T coo = T c0 (A = A(l-d)). (72) 

In the case of (i-wave pairing from Eq.([F(]) we get: 

T coo =T c0 (A = A(l-a d )) (73) 

where 

sin7ro; 

a d = a H (74) 

7T 

defines the "effective" size of flat parts of the Fermi surface for the case of d-wave pairing. 
Thus, for T < T CQO superconducting gap is nonzero for any values of D and drops form the 
value of Ao to Aoo with the growth of D. For T coo < T < T c q the gap is nonzero only 
for D < D max . Appropriate dependences of A on D can easily be obtained by numerical 
solution of Eqs. (|66|) and ( |67| ) . 

In our model of the pseudogap state dielectric gap D is not fixed but random and dis- 
tributed according to (|29|). The above results should be averaged over these fluctuations. It 
is important that within this model we can directly calculate and exactly calculate super- 
conducting gap A, averaged over the fluctuations of D: 

<A>=y o dDV(D)A(D) = — ^ dDDe-^A(D) (75) 



The dependences of A(D) described above lead to the average gap ([75]) being formally 
nonzero up to the temperature T = T c0 , i.e. the temperature of superconducting transition in 
the absence of pseudogap anomalies. At the same time, the temperature of superconducting 
transition T c in a superconductor with the pseudogap is, obviously, smaller than T c0 [|63] . 



This paradoxical behavior of < A > apparently can be explained by D-fluctuation induced 
appearance of local regions in the sample with A ^ (superconducting "drops") in the 
whole temperature region of T c < T < T c0 , with superconducting coherence appearing only 
for T < T c . Of course, the full justification of this qualitative picture can be obtained only 
after the analysis of more realistic model with finite correlation length £ of AFM fluctuations 
23 . At the same time, simplicity of our "toy" model with £ — > oo allows to get an exact 
solution for < A > immediately. 



23 Qualitatively this situation reminds the appearance of inhomogeneous superconducting state, induced by 
strong fluctuations of local electronic density of states, close to the Anderson metal - insulator transition 



[101, 108| 
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To determine superconducting transition temperature T c in a sample as a whole we may 
use the standard mean - field approximation procedure to average over random short - range 
order fluctuations (similar approach is used in the theory of impure superconductors [|101|| ), 
which assumes the self - averaging of superconducting gap over fluctuations of D (i.e., in 
fact, independence of A on fluctuations of D). Then the equations for this mean - field A m f 
have the following form: 

A [ dda-^ dDDe-^-^= 2T +(1-5) , 2T } (76) 







for the case of s-wave pairing, and 



4 r° I 2 roc D , f * ^ yg+ggSl 

- .A, h Jo vry . + D 2 + A 2 mf e 2 ^) 

+ / d0e 2 (0) 



(77) 



for the case of <i-wave pairing. These equations are the limiting case (for £ — > oo) of 
superconducting gap equations analyzed above (|4"3| ) using the recurrence procedure (p9|). 

From (]76f) and (|77|) it is easy to get the appropriate T c - equations, giving the temperature 
at which the homogeneous gap A m f first appears everywhere in the sample. For the case of 
s-wave pairing we have: 

, . \/t 2 +D 2 




In case of d-wave pairing we have just to replace a in ([78]) by "effective" ad defined in ([HI). 
These equations for T c coincide with coincide with those obtained during the microscopic 



derivation of Ginzburg - Landau expansion in this model |)3| and with £ — > oo limit of 
equations obtained using fl3"5|), (|3|). In general we always have T coo < T c < T c0 . 

Temperature dependences of average gap < A > and mean - field gap A m f, obtained 
by numerical solution of our equations for the case of s-wave pairing, are shown in Fig.[Z7| 
24 . We see that A m f goes to zero at T = T c < T cQ , while < A > is nonzero up to 



24 In case of o?-wave pairing temperature dependences of < A > and A m y are qualitatively similar to the 
case of s-wave pairing. 
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FIG. 27: Temperature dependences of superconducting gaps A m j (dots), < A > (full curves) 



and Aq (dashed curves) for the case of s-wave pairing [|102| . 1. - A = 0.4; a = 2/3; u> c /W = 3 



(T c /r c0 = 0.42). 2. - A = 0.4; a = 0.2; u c /W = 1 (T c /T c0 = 0.71). 

T = T c0 and characteristic "tails" in the temperature dependence of < A > in the region of 
T c < T < T c0 , in our opinion, correspond to the abovementioned picture of superconducting 
"drops" appearing in this temperature interval, while superconductivity in a sample as a 
whole is absent. Note that temperature dependences of < A(T) > shown in Fig.^7| are 
quite similar to those observed in underdoped cuprates via ARPES ^9f or specific heat 



measurements and shown in Figs, y, [14], if we assume that T c observed in these samples 
corresponds to our mean - field T c , while "drops" with < A >^ exist in the region of 
T > T c up to T c o, which is much higher than T c . Non self - averaging nature of the gap 
in our model is also reflected in the fact that in the whole temperature region of T < T c q 
gap dispersion < A 2 > — < A > 2 is different from zero, circumstantially confirming the 
proposed qualitative picture of "drops" . Of course, full justification of this picture is possible 
only with the account of finite values of £. Note that the value of T c0 is not very well defined 
from the experimental point of view. Speaking above about superconducting transition 
temperature in the absence of the pseudogap we assumed the value of T c0 of the order of T c 



observed for optimal doping. However, to compare the results shown in Fig.|27] with data of 
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Figs.^, [14| we need the values of T c q much higher than optimal T c . in this sense the described 
picture of superconducting "drops" induced by "dielectric" gap fluctuations becomes not so 
very different from precursor pairing picture of superconducting pseudogap scenario. At the 
same time it is clear that the most part of data shown in Figs. |3|, |n] may be explained by 
the existence of the pseudogap of dielectric nature, not by just discussed non self - averaging 
nature of superconducting A. Note, for example, that superconducting gap obtained from 



tunneling data of Refs. |2T], |22] (appearing on the "background" of much wider pseudogap) 
demonstrates "normal" temperature behavior going to zero at T = T c . 

Despite the fact that superconductivity in a sample as a whole for T c < T < T c0 is absent, 
the existence in this region of nonzero average < A > leads to anomalies in a number of 
observables, e.g. in tunneling density of states or spectral density measured via ARPES 

In particular, the tunneling density of states for the case of d-wave pairing has the fol- 
lowing form [ 102|| : 



A ( E ' 4 2 f °° dDDe-^ { f dxj> , |E| 6{E 2 - A 2 ( J D)e 2 (0) - D 2 ) + 



iVo(O) ttW 2 Jo [Jo Je 2 - D 2 - A 2 (D)e 2 (0) 



+ r /4 d( p |E| Q{E 2 - A 2 {D)e 2 {<P)) 

IE 2 - A 2 (D)e 2 (</>) 



a 



(79) 



Assuming the self - averaging nature of the gap we have A = A m f and independent of 
D. In this case the width of superconducting pseudogap in the density of states is of the 
order of A m / and appropriate contribution disappears for T — > T c , only pseudogap due to 
AFM fluctuations remains. In reality we have in ([79]) A = A(D), which is determined from 

Tunneling density of states for the case of (i-wave pairing is shown in Fig.^. We see 
the significant difference of the exact density of states from that obtained in the mean - 
field approximation, which is due to fluctuations of superconducting gap (superconducting 
"drops") induced by fluctuations of AFM short range order. The exact density of states in 
fact does not "feel" superconducting transition in a sample as a whole at T = T c . Charac- 
teristic scale of superconducting gap (pseudogap) in tunneling density of states becomes of 
the order of A , not A m f, as in the mean - field approximation over fluctuations of short 
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FIG. 28: Density of states in case of <i-wave pairing [102 1 . A = 0.4; a = 2/3; oj c /W = 5 (T c /T c q = 
0.48, T coo /T c0 ~ 10~ 18 ), T/T c0 =: 1. - 0.8; 2. - 0.48; 3. - 0.1. Dotted line shows the mean - 
field density of states N m f(E) for T/T c q = 0.1. Dashed line - density of states with pseudogap for 
T > T c0 . 



- range order. Superconducting contributions to the density of states become observable 
already for T = T c0 > T c . 

These results may, in principle, explain unusually high values of the ratio 2A/T C , observed 
in some tunneling experiments in underdoped cuprates p0| , |109|| , as well as in ARPES ||110 



The presence of superconducting "drops" may also explain unusual diamagnetism of these 



systems sometimes observed above T c [pT 



CONCLUSION. PROBLEMS FOR THE FUTURE. 



Finally let us discuss some general conclusions and remaining problems. Scenario of pseu- 
dogap formation, based on the picture of fluctuations of AFM (SDW,CDW) short - range 
order leads to general qualitative agreement with main experimental facts on pseudogap in 
HTSC - cuprates. Our theoretical analysis was intentionally semi - phenomenological and we 
just "modeled" pseudogap anomalies of electronic spectrum, described by parameters of ef- 
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fective pseudogap width W and correlation length £, which, in principle, may be determined 
from the experiment. Both "hot spots" and "hot patches" models can be solved "nearly 
exactly" [13, [DJ, which allows significant progress in the analysis of pseudogap influence 
upon superconductivity. Let us stress once again that it is not very important whether we 
use AFM - fluctuations model (as most popular in the literature) of the model of short -rang 
order fluctuations of some other nature, leading to partial dielectrization of the spectrum 
(such as CDW or structural). 

Deficiencies of our models are due to simplifying assumptions necessary to obtain these 
"nearly exact" solutions. The main of these are the static approximation and assumption 
of Gaussian nature of short - range order fluctuations. The account of dynamic nature of 
fluctuations is absolutely necessary for low enough temperatures, particularly in supercon- 



ducting phase where pairing itself may be due dynamics of these fluctuations |58|, p9|| . We 
feel, however, that simplified analysis given above may in fact describe the most impor- 
tant effects of pseudogap renormalization of electronic spectrum (pseudogap formation on 
"hot" parts of the Fermi surface) and its influence upon superconductivity. The account 
of dynamics of spin fluctuations will inevitably lead us outside the simple phenomenology 
of BCS model. Assumption of Gaussian statistics of fluctuations also may be justified only 
for high enough temperatures and not very close to the line of antiferromagnetic instability. 
However, even this simplest approximation allows qualitative description of almost all basic 
anomalies due to pseudogap formation. Without this assumption the simple and physically 
appealing structure of our equations is destroyed, but only this structure allows us to analyze 
superconducting phase in a relatively simple manner. 

The extension of our semi - phenomenological approach can be apparently done within 
the framework of full microscopic analysis of the problem, e.g. on the basis of the Hubbard 



model. Such attempts were made in a number of papers, e.g. in Refs. 11, 12, 1 3 1 cited 



above, as well as in Refs. ||112| , |113|| . Though many qualitative conclusions of these papers 
coincide or are close to those discussed above, the analysis is in most cases restricted to first 
diagrams of perturbation series ||114| , |115|| and some kind of self - consistent procedure on 
this basis. In this case there is no big problem also to take into account also the dynamics 
of spin fluctuations. However, any generalization to higher order contributions is difficult 
enough. 

Microscopic justification of the existence of the wide region of the phase diagram with 
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well developed fluctuations of AFM (SDW) short - range order is, of course, of principal 
importance. It should be noted that this region of "critical" fluctuations may be wide enough 
just due to the low dimensionality (quasi two - dimensional nature) of HTSC - systems. For 
example, the region of superconducting critical fluctuations in HTSC - cuprates may reach 



the width of some tens of degrees K [|116 . As AFM fluctuations are characterized by an 
energy scale which is an order of magnitude higher, the existence of the "critical" region with 
the width of some hundreds degrees does not seem improbable. However, the microscopic 
justification of such picture is still absent. 

Probably the main qualitative conclusion from our theoretical picture concerns the old 
discussion about Fermi - liquid or non Fermi - liquid nature of electronic spectrum in HTSC 
- cuprates. According to our simple models the electronic spectrum (spectral density) is 
Fermi - liquid like only on "cold" parts of the Fermi surface (in the vicinity of Brillouin zone 
diagonals), while on "hot" parts the spectral density is in general non Fermi - liquid like due 



to strong scattering on AFM fluctuations |T4|, [15| . The presence or absence of Fermi - liquid 
behavior essentially depends on the value of correlation length of short - range order £ (Cf. 



Refs.[81, 82|). In a recent paper | 117 | this picture was criticized on the basis of successful 
enough fit of experimental ARPES data of Refs.[|51], |118| , |119|| for optimally doped Bi — 2212 
to dependences of the "marginal" Fermi - liquid theory. Particularly, in Ref . ||1 1 7|| it was 
claimed that Fermi - liquid like behavior of the quasiparticle damping is absent everywhere 
on the Fermi surface and there is no significant momentum dependence of this damping. 



In fact the authors of Ref. [|117|1 had to introduce the significant momentum dependence of 
static scattering rate attributed to anisotropic impurity scattering, which was necessary to fit 
data, even for optimally doped system. Linear dependence of damping rate on the energy of 
quasiparticles, which is observed in HTSC - cuprates (and postulated in "marginal" Fermi 
- liquid theory), can not disprove the usual Fermi - liquid theory, because the standard 
quadratic dependence of da [21 mping may be observable only in some very narrow energy 
interval close to Fermi level, which is apparently smaller than typical precision of ARPES 
experiments. Rather detailed discussion of this point of view can be found in recent reviews 
pl| , |93|j . It must be stressed that the existence of anomalous scattering by AFM fluctuations 
(with wave vectors of the order of the vector of antiferromagnetism) is beyond any doubts 
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2r, 



Our theoretical discussion also used rather traditional scheme of averaging over the ran- 
dom field of AFM (SDW,CDW) fluctuations, which assumes the average spatial homogeneity 
of the system. At the same time a number of experimental data and theoretical models stress 
the possible phase separation in some HTSC - cuprates (especially in underdoped region) 
| |121| , |122|| . This phase separation takes place on microscopic scales and the system consists of 
"metallic" (superconducting) and "dielectric" (antiferromagnetic) "stripes" with the width 
of the order of few interatomic distances. Naturally this situation can hardly be described 
using standard methods used above, while in the opinion of many authors these "stripes" 
may be of primary importance for the physics of HTSC. We shall only note, that the quali- 
tative picture of the random field of AFM (SDW,CDW) fluctuations introduced above also 
implicitly assumes the appearance of effectively antiferromagnetic regions of the size of the 
order of ~ £ intermixed with the regions of the same size, where AFM order is destroyed. In 
this sense the picture used above seems to be not drastically different from the concept of 
phase separation ("stripes"), being different mainly by the smoothness of borders between 
the "metallic" and "non - metallic" regions. Note also that the formalism of the number of 
models of phase separation is based upon the idea of the closeness of the system to some 
kind of CDW instability p3, IT23, |124|| and explicitly introduces scattering by fluctuations 



with special wave vectors. 

It is well known that a lot of anomalies is observed in HTSC - systems under strong 
structural disordering [|101|| . The question of the importance of structural disordering in the 



pseudogap state is almost not investigated. We may mention only few papers, where the 



effects of random impurities were studied ||125| , 126 . From these preliminary results it may 



be concluded that controlled disordering may be rather informative method of studies of the 
pseudogap state. Theoretically these problems are not studied at all. 

The studies of pseudogap anomalies in HTSC - systems can lead also to rather important 
conclusions with respect to possible applications. We have seen that pseudogap formation 
in some sense suppresses superconductivity and, in fact, such important characteristics of 



25 Note an interesting recent paper [ 120 1, where rather good description of unusual transport properties 
of HTSC - oxides was obtained within the model of simplest possible Boltzmann like scattering by spin 
fluctuations, assuming that the contribution of "hot" parts of the Fermi surface to electron transport is 
just excluded. 
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superconductors as critical current and critical magnetic fields are maximal not at optimal 
doping po w 0.15 — 0.17, corresponding to maximal T c , but in the vicinity of the "critical" 
doping p c pa 0.19 ||127|| . This fact may be of some significance for optimization of technical 
superconductors made of HTSC - cuprates. 
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